Two-Point Entanglement Near a Quantum Phase Transition 
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In this work, we study the two-point entanglement S{i,j), which measures the entanglement 
between two separated degrees of freedom (ij) and the rest of system, near a quantum phase 
transition. Away from the critical point, S{i,j) saturates with a characteristic length scale ^e, as 
the distance \i — j\ increases. The entanglement length agrees with the correlation length. The 
universality and finite size scaling of entanglement are demonstrated in a class of exactly solvable 
one dimensional spin model. By connecting the two-point entanglement to correlation functions in 
the long range limit, we argue that the prediction power of a two-point entanglement is universal 
as long as the two involved points are separated far enough. 
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Quantum phase transition happens at zero tempera- 
ture when physical parameters are changed^. Just like 
the classical critical phenomena^, a diverging length scale 
dominates the physics near a quantum phase transition. 
Quantum criticality also exhibits scaling law and uni- 
versality class. Meanwhile, quantum system contains a 
new physics, entanglemenlj^. Entanglement is a unique 
property of quantum system and measures quantum cor- 
relations. It is thus natural to ask if entanglement pos- 
sesses the dominating correlations near a quantum phase 
transition. Tremendous interests and work have been in- 
vested along this line. In particular, entanglement of for- 
mation between two separate spins after tracing out the 
rest has been studied widely^'^i^'ii^'^. It is shown that 
this quantity shows pronounced behaviors such as scal- 
ing law and universality class^i^. However, this quan- 
tity vanishes as two spins are separated by several lat- 
tice constants. Thus, it does not contain the dominat- 
ing long-range feature between two separated spins pos- 
sessed by the system. The entanglement between a block 
of spins and the rest of system is also studied in vari- 
ous modelsiSiiii'^ ? ? ? However, there is no direct re- 
lation between it and the long range correlation of the 
system either—. Recently, Verstraete et al^^ studied the 
average entanglement that can be localized between two 
separated spins by performing local measurements on 
the other individual spins. The typical length scale at 
which the localizable entanglement decays is defined as 
the entanglement length. All classical correlation func- 
tions provide lower bounds for localizable entanglement. 
However, the entanglement length defined for localizable 
entanglement can be infinite even for a gapped system^^. 
Therefore, an entanglement measure describing the cor- 
rect long range physics near a continuous phase transition 
has yet to be studied. It is the purpose of this work to 
identify and study such a measure of entanglement near 
a quantum phase transition. 

We start by noticing that two-point correlation func- 
tions can be essential in the study of quantum phase tran- 
sitions. It is also convenient to organize correlation func- 
tions into a reduced density matrix. We will thus focus 
on the possible entanglements encoded in a reduced den- 
sity matrix. The reduced density matrix of two separated 



points i and j contains at least three types of entangle- 
ment. The first type is the entanglement between i (or 
j) and the rest of system. It can be measured by the 
Von Neumann entropy of the reduced density matrix of 
single point i, which in turn can be obtained from the 
two-point reduced density matrix by partial tracing out 
j. This quantity has been studied near quantum phase 
transition^ and does not contain an explicit length scale. 
The second type is the entanglement between i and j and 
can be measured by concurrence for spin-1/2 systems. It 
shows beautiful scaling law but does not contain a di- 
verging entanglement length scale near a quantum phase 
transition^. In this work, we shall focus on the behav- 
ior of the third type of entanglement, the entanglement 
between (ij) and the rest of system. At zero temper- 
ature, the whole system is described by a pure state. 
It is thus justified to use the Von Neumann entropy of 
the two-point reduced density matrix to measure the en- 
tanglement between the two points (ij) and the rest of 
system. We shall call this type of entanglement as two- 
point entanglement. It is shown that two-point entan- 
glement is dominated by a diverging length scale in the 
vicinity of a quantum critical point. We show the en- 
tanglement length diverges as approaching a quantum 
phase transition and agrees with the correlation length 
in an exactly solvable spin model. We also demonstrate 
that two-point entanglement exhibits beautiful scaling 
law and universality2iii. The non-analyticity of the two- 
point entanglement for two nearest-neighboring points 
has been successfully utilized to efficiently identify pos- 
sible phase transitionsi^iiiii^ii^iil. Because of the nice 
scaling properties, phase transition point can be deter- 
mined from small systems with considerable accuracy 
without pre-assumed order parameters^ii. However, stud- 
ies along this line are so far limited to specific models and 
it is not clear how universal its prediction power of a pos- 
sible phase transition is. By connecting the two-point en- 
tanglement to correlation functions, we shall argue that 
the prediction power of a two-point entanglement is uni- 
versal as long as the two involved points are separated 
far enough. The tow-point entanglement has great ad- 
vantages in predicting possible phase transitions. Al- 
though it does not provide extra information beyond cor- 
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relation functions, two-point entanglement automatically 
picks out the most dominant non-analyticity of all possi- 
ble correlation functions. Therefore, one single quantity 
is enough to distinguish phase transitions between differ- 
ent orders. Pre-assumed order parameters are thus not 
needed. This is particularly convenient when studying 
systems with complicated phase diagram, for instance, 
strongly correlated systems such as Hubbard modeii^. 
We also want to point out that two-point entanglement 
is different from the entanglement entropy of block spins 
studied widely in literature. The reduced density matrix 
of two spins at i and j contains only a subset of the infor- 
mation encoded in the reduced density matrix of a block 
of spins with size of the block L > \i — j\. Because of 
the non-additive nature of the entanglement entropy, the 
behavior of two-pint entanglement studied in this work 
can be very different from the one of block spins. The 
entanglement entropy of a block of spins contains more 
information than our two-point entanglement does. How- 
ever, it is one of our main results that the two-point en- 
tanglement bears enough information to efficiently iden- 
tify possible quantum phase transitions. The two-point 
entanglement is also closely related to the mutual infor- 
mation recently studied by Gu et al'^ . In a translation- 
ally invariant system, two-point entanglemtn and mutual 
information share the same range-dependence and thus 
bare the same characteristic length scale. However, two- 
point entanglement has an advantage that it is a better 
measurement of quantum information in the sense that 
the mutual information shared by two systems may be 
larger than the entropy in either part. 

As an example, let us study the two-point entangle- 
ment in an spin- 1/2 ferromagnetic chain with an ex- 
change J and a transverse field h, 
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where (T^,cr^,cr^ are Pauli matrices, N is the number of 
sites and 7 is the anisotropy in the xy plane. This model 
belongs to the Ising universality class when < 7 < 1 
and N = 00. It can be solved exactly and the correlation 
functions are given in Reffi^ii^i^i^. 

The two-point entanglement S(i,j) is defined as the 
entropy entanglement between two separate spins, i and 
j, and the remaining spins in the system. From the full 
density matrix of the whole system, we can obtain the 
reduced density matrix p for i and j by partial tracing 
out the rest spins. Based on the symmetry of the Hamil- 
tonian, we know that p is real and symmetric^. The 
nonzero elements are pu, P22 — P33, P44, P14 and p23 if we 
choose the basis |1) = | TT), |2) = | Ti), |3) = 1 it), |4) = 
I II). These elements can be related to the correlation 
functions p"{i,j) = (afcr") — {cr'j'){(Jj) for a = x,y,z 
and magnetization M{i) — {pf}^. We can then define 
the two-point entanglement as the entropy of p 




FIG. 1: S(n) as a function of distance n for different values of 
A in the Ising limit where 7 = 1. As A starts from small value 
and approaches the critical value 1, it takes longer and longer 
distance for 5(n) to saturate. A length scale, entanglement 
length ^_E, can be defined to measure how fast the two-point 
entropy saturates. 



S'(«, j) measures the entanglement between two separated 
spins, i and j, and the rest spins. Since the system is 
translationally invariant, we know and S only depen- 
dent on the distance n = \i—j \ and M is site-independent. 

Using the results of Ref ]^^i^^'^°i^^ , we can calculate the 
reduced density matrix p and thus S exactly. Let us in- 
troduce the dimensionless tuning parameter A = J/h. 
There is a quantum phase transition at Aj, = 1 for an 
infinite spin chain. For A > 1, the system is ferromagnet- 
ically ordered and kinks start to condense at A = 1 which 
drives the system into a quantum paramagnetic state in 
A < 1 region. 

In FiglU we show S{n) as a function of distance n for 
different values of A in the Ising limit where 7=1. S{n) 
increases as n increases. This is opposite to the decaying 
of correlation functions. We want to emphasize that S 
measures the entanglement between two spins (ij) and 
the rest spins while correlation function measures the 
correlation between i and j. Crudely speaking, S can 
be taken as an accumulation of the correlation attached 
to i or/and j. Consequently, the two-point entanglement 
S(n) is expected to increase and reach a constant as n 
increases. FigH] reveals a more important message. As 
A starts from a small value and approaches the critical 
value Ac = 1, it takes longer and longer distance for S(n) 
to reach the long-range limit S{oo). One can thus define 
a length scale which measures how fast the two-point en- 
tanglement S saturates. This can be shown by examining 
the long range behavior of S{n) for A = Xc{N) — 6X near 
Xc{N). In the thermodynamic limit N 00, the long 
range behavior of correlation functions has been stud- 
ied in Rcf.~iS2ii^. From Ref.^^, the asymptotic forms of 
correlation functions for A < 1 and large n are 



SiiJ) = -TrpXogp. 



(2) 
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(3) 
(4) 
(5) 
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with 

In the hmit of 7 ^ 1, we have A2 A. Near the critical 
point Ac = 1, A2 « 1 + (A — l)/-f to the first order of 
A — 1. Since and vanish faster than as n goes to 
infinity, it is safe to set them to zero for large enough n. 
The eigenvalues of the reduced density matrix p are thus 

= (l-Af2±p-(n))/4 (7) 

It follows that 

Sin) - S'(oo) oc oc n-^e-^^^^ (9) 

where the entanglement length is given by 

(10) 

Eq.® directly relates the entanglement to the correla- 
tion function. The expression Eq.Q is essentially a Tay- 
lor expansion of S{n) in terms of the correlation functions 
and thus holds as long as the correlation functions are 
small, i.e., when n is longer than the correlation length. 
In Figl2] (a), we compare the exact results of S{n) and 
the scaling relation given by Eq.® for A = 0.9. The plot 
of ln(S'(n) — S{oo)) vs n for A — 0.9 confirms the entan- 
glement length given by Ea. pO| directly. It is interest- 
ing to notice that is directly related to the correlation 
length. The critical exponent v = lis also independent of 
the anisotropy 7, in agreement with the universality hy- 
pothesis. Similarly, from the asymptotic forms for A > 1 
given in Fig.l of Ref.^^, we again obtain oc (A — 1)^^ 
with critical exponent v = \. Eg. pH)) suggests that the 
entanglement length diverges at the critical point A = 1 
and the two-point entanglement should approach S{oo) 
through a power law. We shall now demonstrate that 
this is indeed the case. For A = 1 and large n, we have^^ 

p^{n) oc 71-1/4 

py{n) cx (12) 
cx (13) 

The two-point entanglement can be shown to follow 

5(n) - 5(00) oc (14) 

where the exponent —1/2 is related to the anomalous 
dimension of . In Fig[2] (b), we compare the exact 
results of S{n) and the scahng relation given by Eq.® 
for A = 1. The plot of ln(S'(n)-5(oo)) vs Inn for A = 0.9 
confirms the exponent —1/2 in the scaling law Eq. (fT4l) 
directly. 

We thus arrive at the main result of this work. The 
two-point entanglement S{n) contains a length scale, the 
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FIG. 2: (a) The comparison between exact results of S{n) 
for A = 0.9 and Eq.® and Eq.dlOj. The ln(S'(n) - S{oo)) 
vs n plot directly confirms the entanglement length defined 
in Eq. (|10|) . Inset is the plot of S{n) as a function of n. (b) 
The comparison of S{n) for A = 1 and Eq. (fT4)) . The plot of 
ln(5'(n) — S{oo)) vs Inn for A = 0.9 confirms the exponent 
— 1/2 in the scaling law Eq. (|14|l directly. Inset is the plot of 
S{n) vs n. 



entanglement length that agrees with the correlation 
length. For < 7 < 1, S{n) reaches 5(00) exponentially 
over the length scale S^e- £.e diverges as A approaches 
the critical value Ac = 1 with an exponent v — I. At 
A = 1, S{n) follows a power law given by Eq. p^ with an 
exponent —1/2. Therefore, the two-point entanglement 
bears the same dominating long range physics about the 
phase transition as correlation functions do. This is not 
surprising since the reduced density matrix p and con- 
sequently the two-point entanglement S{n) are smooth 
functions of correlation functions. We also notice that 
the entanglement length S.E defined in this work is closely 
related to the one defined in Re£— , where it is defined 
as the measure of how fast the localizable entanglement 
between two separated spins decays as the distance in- 
creases. In a very crude model, as we mentioned previ- 
ously, the two-point entanglement S{i,j) can be taken as 
an accumulation of entanglement. Only those spins that 
are entangled with one or both of i and j can contribute. 
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n=2, N=1001 



FIG. 3: Top: S{n); bottom: dS{n)/d\ as functions of A for 
different n. For finite n, tfie maximum of S does not corre- 
spond to a phase transition wfiile sharp dip or peak in dS/dX 
does. When n oo the maximum of S shifts to the critical 
point and displays a A-like anomaly. The anisotropy 7 is set 
to 1 in this plot. 



Therefore, through the study of the distance dependence 
of the two-point entanglement, we learn how the entan- 
glement is shared out among spins. If the localizable en- 
tanglement attached to either spin decays exponentially 
with a length scale ^, one then expects S(n) increases 
as n increases but still shorter than ^. After n >> ^, 
S(n) stays at a constant as further separation no longer 
introduces extra localizable entanglement. However, the 
contribution of localizable entanglement to the two-point 
entanglement is in general not greater than the optimized 
one defined in Refi^^. Thus, the two definitions of en- 
tanglement length can be very different in some cases. 
For instance, it has been shown that the decay length 
of localizable entanglement can be infinite in a gapped 
systemSS. In contrast, the entanglement length defined 
in this work is expected to follow the correlation length 
which is finite' . 

To further explore the connection between two-point 
entanglement and quantum phase transition, we plot 
S(n) and dS{n)/d\ as functions of A for different n in 
FigO As shown previouslySiil, the maximum of S{n) 
for finite n does not signal a possible phase transition. 
The sharp peak (n = 1) or dip (n > 4) in dS/\ does 
correspond to a phase transition. More interestingly, 
when n increases, the maximum starts to shift towards 
the critical point and develops a A-like cusp at the crit- 
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FIG. 4: (a) The second derivative of <S'(2) with respect to A. 
This plot confirms the kink of 5(2) near A = 1, although the 
anomaly in this case is much weaker than ^(I). Left inset 
is S and right inset is the first derivative, (b) The second 
derivative of S(n) for n = 1,3,5,17. The non-analyticity 
near A = 1 is evident in this plot. 



ical point when n becomes infinity. Another interest- 
ing point is that as n increases, the anomaly becomes 
stronger and stronger as suggested by the larger and 
larger \ds / d\\\^x^. This is another evidence that the 
two-point entanglement catches the important long range 
physics near a phase transition. The short range two- 
point entanglements also exhibit anomaly near the criti- 
cal point. In particular, the curve with n = 1 displays a 
peak at critical point while long range ones have a dip. 
For n = 2, 3 the anomalies near the critical point are 
very weak showing crossover from strong peak of n = 1 
to strong dip for n > A. We also want to emphasize that 
5(2) is non-analyticity at A = 1 although it is not as ap- 
parent as S{1) in the first derivative plot. To show this, 
we plot the second derivative in Fig|4l^a) , which confirms 
a kink near A = 1. In FiglUb), we also plot the sec- 
ond derivative for different n. It is straightforward to see 
the non-analyticity near A = 1 in this plot. With this 
plot, the difficulty of identifying phase transition from 
n — 2, 3 is circumvented. The entanglement <S'(1) in- 
volves the correlation function {(^fo'i^i), which is coin- 
cidently the expectation value of the coupling term in 
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FIG. 5: Finite size scaling of dS{7)/dX with exponent f — 1. 
The calculation was carried out for systems of difference sizes, 
varying from A'' = 61 to A'^ = 601. Left inset: The first 
derivative at Xc{N) diverges logarithmically; Right inset: The 
position of the dip moves toward the true critical point Ac = 1 
as the system size A'' increases. Again, 7 = 1 in this plot. 



the Hamiltonian^^. This imphes that the peak anomaly 
of dS{l)/dX is dominated by short range interactions of 
the Hamiltonian. One thus might expect that its pre- 
dicting power of a possible phase transition is limited 
and depends on the details of interactions. On the other 
hand, we have argued that the non-analyticity in long 
range two-point entanglements is dominated by a diverg- 
ing length scale £,e and is expected to be universal. 

Typical finite size scaling for dS{n)/dX is shown in 
FigEl The data of dS{7)/dX collapse onto a universal 
line for different sample sizes N and different coupling 
strength A. The plot contains data from system size vary- 
ing from = 60 up to = 701. The critical exponent 
extracted from the scaling is again v = I, in agreement 
with our previous discussion that the two-point entangle- 
ment S{n) contains the correct physics dominated by a 
diverging length scale near the critical point. In the left 
inset of Figini we also plot the dip value dS{7)/dXc which 
diverges logarithmically. The position of the dip moves 
toward the true critical point Ac = 1 as the system size N 
increases. It has been shown recently that dS{l)/dX also 
bears beautiful finite size scaling with exponent ly = 1^^. 
We have also verified that similar scaling law holds with 



critical exponent v = I when < 7 < 1 is considered. 
As an example, we plot the finite size scaling of dS{3) / dX 
for 7 = 0.5 in FiglB] 

In summary, we study the two-point entanglement 
S{n) which measures the entanglement between two sep- 
arated degrees of freedom and the rest ones of the system. 
We introduce the entanglement length that measures 
how fast the entanglement saturates as the distance in- 
creases. Two-point entanglement automatically picks out 
the most dominant non-analyticity of all possible two- 
point correlation functions. Consequently, its entangle- 
ment length agrees with the correlation length of the 
most dominant correlation functions. The universality 
and finite size scaling of entanglement are also studied in 
a class of exactly solvable spin model. 




FIG. 6: Finite size scaling for dS{3)/dX with 7 = 0.5. The 
size A'^ varies from 61 up to 1001. Inset plots dS{n)/dX for 
different n with A^ = 401. 
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